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Transverse- field muon spin rotation measurements of overdoped La2-xSrxCu04 reveal a large 
broadening of the local magnetic field distribution in response to applied field, persisting to high 
temperatures. The field-response is approximately Curie- Weiss like in temperature and is largest for 
the highest doping investigated. Such behaviour is contrary to the canonical Fermi- liquid picture 
commonly associated with the overdoped cuprates and implies extensive heterogeneity in this region 
of the phase diagram. A possible explanation for the result lies in regions of staggered magnetization 
about dopant cations, analogous to what is argued to exist in underdoped systems. 

PACS numbers: 74.25.Ha; 74.62.Dh; 74.72.Dn; 76.75.+i 



I. INTRODUCTION 

Of the many recurring themes in cuprate research, 
one of the most persistent is atomic-scale heterogene- 
ity. In hole-doped cuprates in particular, local probes 
have provided numerous examples where superconduct- 
ing, electronic and magnetic properties vary on a length- 
scale of tens of angstroms. Scanning tunneling spec- 
troscopy (STS) regularly reports variations in the su- 
perconducting gap of Bi2Sr2CaCu208+5 (BSCCO) which 
see m to cor relate well with the location of dopant oxygen 
^Qj^^TmnEEl^ Similar measurements on La2-xSrxCu04 
(LSCO) reveal a spatial variation in the local density 
of states^. A number of NMR and nuclear quadrupole 
resonance (NQR) studies have reported inequivalent 
sites an d doping he terogeneity in La-based cuprate 
system&^EEn™™^^'^ 

Magnetic properties of underdoped systems are domi- 
nated by incommensurate order and fluctuations, com- 
monly associated with one-dimensional variations in 
charge and spin density^^. Zero- field muon spin rota- 
tion (ZF-/iSR) measurements of Lai.88Sro.i2Cu04 and 
La2CuO4.11 suggest that such heterogeneous 'stripe or- 
der' can itself be fragmented, and that long-range 
order in these systems results from percolation be- 
tween overlapping 15- 30 A island^^. Separate studies of 
Lai.85-yEuySro.i5Cu04"^^ and La2-xSrxCu04+y^ report 
similar phase separation over a wide range of doping, 
with superconducting and magnetic volumes mutually 
exclusive and oppositely correlated. 

Upon doping Zn^+ ions onto the planar Cu^+ sites of 
YBa2Cu306+5 (YBCO), a number of NMR studies re- 
port the existence of staggered spins immediately about 



the substitution ^^t^^^H^^I^m, A local 
suppression of superconductivity in the same regions of 
the sample has been inferred from first /iSR^^ and later 
STS^^ studies, though specific interpretatio n of both 
NMR and /iSR results is still a matter of debat^^^. Re- 
cent theoretical work suggests that similar antiferromag- 
netic islands are expected in underdoped LSCO around 
dopant Sr^+ ions due to correlation effects, and that this 
may explain several properties of the La-based cuprate 
superconductors^^ . 

Whereas heterogeneity is well-established in the un- 
derdoped cuprates, the situation is much less clear on 
the overdoped side of the phase diagram. It is widely 
believed that the overdoped materials behave like canon- 
ical Fermi-liquids, and this belief is supported by select 
transport measurement Effects from Zn^+ 
impurities in YBCO are also less pronounced beyond op- 
timal doping^^ 22 | 23 | 25 |^ QYen glassy order disappears 
around x=0.19 in LSCO^. On the other hand, bulk sus- 
ceptibility studies of both LSCO and Tl2Ba2Cu06+5 (Tl- 
1201) observe a Curie component emerging with progres- 
sive overdoping^^ , in stark contrast to the tem- 
perature independent Pauli susceptibility expected for 
a Fermi liquid. Moreover, STS on heavily overdoped 
Bi2-xPbxSr2CuOy reports nanoscale variations in the 
superconducting gap^, unpaired carriers have been ob- 
served at low temperatures in overdoped LSC033, and 
real space phase separation has been argued indirectly 
from /iSR^^, optical conduct ivit}^^^^, magnet iz at ioiPI' 
and neutron scattering^^ measurements. 

In recent years, a handful of transverse-field (TF) 
/iSR studie^^^^^ have reported an heterogeneous mag- 
netic field response in single crystals of underdoped 



2 



LgCOPEm and YBCCf^, in the form of a symmet- 
ric broadening of the local field distribution in response 
to applied field. First reported by Savici et al. in crys- 
tals containing fragmented stripe ordei^^, the effect was 
originally associated with the staggered susceptibility of 
the heterogeneous spin system. Data from Ishida et 
al]^ at a lower doping were interpreted in a similar way. 
Most recently, however, Sonier et al.^^ demonstrated that 
the size of the field-induced broadening closely correlates 
with superconducting Tc in underdoped YBCO, and in- 
ferred local patches of superconductivity at temperatures 
well-above the bulk transition. 

In this paper, we address the magnetic field response 
of overdoped LSCO with TF-/iSR. Preliminary work by 
us suggests a significant response to applied field at very 
high dopings^^. Here, we present data on four crystals 
ranging from optimal to heavily overdoped in a variety 
of applied fields and temperatures, comparing directly to 
previous work. We show that in every crystal investi- 
gated there exists a spatially heterogeneous response to 
applied field which increases with decreasing temperature 
in a Curie- Weiss-like fashion. The dominant trend is a 
linear increase of field response with x through the origin, 
with a notable change in temperature dependence above 
optimal doping. In Lai.7Sro.3Cu04, we further show that 
this response is highly anisotropic, being larger when the 
field is applied perpendicular to the copper-oxide planes. 
Several possible pictures are put forth, and the ability of 
each to explain our data is discussed. 



II. EXPERIMENTAL DETAILS AND RESULTS 
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FIG. 1: (a)-(e)RRF asymmetries resulting from TF — fiSR 
measurements of single crystal Lai.7Sro.3Cu04 in an applied 
field of /ioH=3T at several temperatures. Field was applied 
perpendicular to the copper-oxide planes. (f)FFTs of three 
time spectra. The vertical dashed line shows the location of 
the applied field. 



Four crystals of LSCO were grown via the traveling 
solvent floating zone method. Crystals with x=0.15 and 
x=0.19 were grown at AIST Tsukuba, and crystals with 
x=0.25 and x=0.30 were grown at the University of 
Toronto. Samples with x=0.19, x=0.25 and x=0.3 were 
cut into thin plates with the c-axis perpendicular to the 
large face. An additional crystal with x=0.30 was cut 
with the c-axis in the plane of the plate to test for signal 
anisotropy. The crystal with x=0.15 was a half cylin- 
der with approximate radius of ~5mm. Crystals for the 
larger two dopings are from the same batch as those stud- 
ied earlier by Wakimoto et al. with neutron scattering^^. 
SQUID magnetization measurements on these two sam- 
ples confirmed results by Wakimoto on related powder 
sampleJ^, where a Curie component grows with over- 
doping by an amount equivalent to 0.5//^ per additional 
hole. 

TF-ziSR^*" data were collected using the HI-Time spec- 
trometer on the M-15 beamline of TRIUMF in Vancou- 
ver, Canada. Fields ranging from 0.2T to 7T were applied 
perpendicular to the large face of the crystals, equiva- 
lent to the crystalline c-axis for all but one sample. Ap- 
plied fields were measured independently by simultane- 
ously observing the precession of muons landing in an 
adjacent crystal of Ca(C03). 



Figures [TJa)-(e) show sample /iSR spectra from 
Lai.7Sro.3Cu04 in an applied field of /ioH=3T ||c at five 
temperatures between 100 and 2K. A dramatic relaxation 
of our TF-/iSR asymmetry is evident, with the rate of 
relaxation increasing more rapidly as the temperature is 
lowered. This is in contrast to separate measurements of 
the same crystal in zero fielcP^, which shows no apprecia- 
ble temperature dependence down to lOK. In the absence 
of magnetic fiuctuations, the TF-/iSR asymmetry func- 
tion is proportional to the cosine Fourier transform of the 
local distribution of fields at the muon sit^^, and thus 
the relaxation rates of the time spectra are a measure 
of the width of the local field distribution at the muon 
site. This is demonstrated explicitly in Fig. [TJf), which 
shows fast Fourier transforms of three time spectra. Of 
particular note is lack of net shift of the average field 
from the applied field (dotted line), despite a significant 
increase in the width of the distribution. Similar spec- 
tra were seen for every crystal investigated and for every 
(non-zero) applied field. 

In order to quantify the magnetic response, we fit all 
spectra to the form: 

A{t) = Ao cosico^t + (^)e-^'''e-^^ (1) 
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where c<;^ is the mean precession frequency of the muons 
(oc average local field), and A is the rate of field-induced 
relaxation. The Gaussian term in Eq. ([T]) is present to 
account for finite field-width due to randomly oriented 
nuclear dipole moments and applied field heterogeneity. 
It is temperature independent and was held constant to a 
value chosen by separate calibration measurements. The 
exponential form for the field response is equivalent to 
a Lor entzian- like broadening of the field distribution and 
was chosen to allow direct c omparison of the current work 
to previous studiej ^ * ^^ * ^^[ Careful examination of the 
/iSR lineshape reveals that spectra are better described 
by a root exponential function at low temperatures, but 
all spectra were well-described by Eq. ([T]) except when 
x=0.30 and T<10K. 
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FIG. 2: (Left) The field dependence of the /iSR linewidth for 
each of the four crystals, extracted from Eq. [l] (Right) The 
slopes of the ^ curves as a function of temperature. Solid 
lines in both panels represent fits described in the text. 

Figure [2] shows the extracted values of A for each of the 



four crystals as a function of applied field and tempera- 
ture. Relaxation rates are largest for our highest doping, 
and generally get smaller with decreasing doping- though 
the x=0.15 sample shows a slightly larger effect than 
the x=0.19 sample. At all temperatures, the rates are 
roughly linear in field and extrapolate to small positive 
values, revealing a slight overestimation of <j in Eq. 
Slopes of the A vs. /ioH curves increase gradually as the 
samples were cooled, with an approximate Curie- Weiss- 
like temperature dependence. 

The temperature dependence for each crystal was fit 
to the form: 



C{x) 



e{x) 



B{x), 



(2) 



where B(x) was included to account for possible tem- 
perature independent sources of relaxation, such as 
Pauli susceptibility from a heterogeneous metallic state, 
Van-Vleck susceptibility from randomly placed dopant 
cations, and poor estimation of the applied field hetero- 
geneity in Eq. [T] Curves of best fit are shown as solid 
lines in Fig. [2] Though Eq. [2] provides a useful character- 
ization of the current data, it is important to note that 
other forms for the temperature dependence describe the 
data equally well. 

The doping dependence of the heterogeneous field re- 
sponse is summarized in Fig. [Sj In addition to the data 
presented above, we have also included data previously 
published in Ref. l53l (x=Q.12). and data extracted from 
Refs. ,54 (x=0.07) and ^ (x=0.145, 0.166 and 0.176). 
The topmost panel shows the increase in ^ ^ at T=50K 
over that at T=100K. The current data set is consis- 
tent with the work of Sonier et al. over the same dop- 
ing range. The overall evolution of the field response in 
the underdoped materials is also qualitatively consistent 
with the trend reported by Sonier in YBCO, where A 
varies in a way that reflects superconducting Tc. Above 
optimal doping however, the system breaks from this 
trend, and the relaxation increases rapidly with Sr^+ 
content up to the largest doping explored. The mid- 
dle and lower panels of Fig. [3] show the respective val- 
ues of C{x) and Q{x) extracted from fits to Eq. [2] Fits 
of data from Refs. [54] and |55| were performed assuming 
a linear field dependence through the origin, which ne- 
glects field-independent sources of relaxation. Though 
there is some scatter, the dominant trend for C{x) is 
a near linear increase with x, passing through the ori- 
gin at zero doping. Fit values for 0(x), large and pos- 
itive in the underdoped materials, go negative at slight 
overdoping and tend towards zero as x get large. We 
note that this critical doping has in the past been asso- 
ciated with the disappearance of glassy order in ZF-/iSR 
measurements^^ and the emergence of Curie term in the 
bulk magnet iz at ion^^^^^J^^^. Specific interpretation of 
these trends will be discussed below. 

Additional measurements were made with a separate 
crystal of Lai.7Sro.3Cu04 with H||6. The main result 
is demonstrated in Fig. [4j The left panel compares the 



4 




T — ' — I — ■ — I — ' — I — ' — I — ' — r 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 



X 



FIG. 3: The doping evolution of the heterogeneous field- 
response. The top panel shows the increase of the absolute 
value of ^# at T=50K over that at T=100K. The follow- 
ing two panels show values of C and G, extracted from fits to 
Eq. |2] The line in the centre panel is a guide to the eye. 



spectra at T=2K, with a field of /ioH=3T applied par- 
allel and perpendicular to the copper-oxide planes. The 
right panel shows the result of fitting such spectra with 
Eq. [T]at several temperatures. A large signal anisotropy 
is immediately apparent, with the heterogeneous field re- 
sponse an order of magnitude larger when the field is 
applied perpendicular to the planes. This anisotropy 
is significantly lar ger than that observed by Savici in 
Lai.875Bao.i25CuOj^. 



FIG. 4: (a)Raw asymmetry from TF — /jSR measurements of 
Lai.7Sro.3Cu04 at T=2K and /ioH=3T, applied perpendicu- 
lar and parallel to the copper-oxide planes. (b)The relaxation 
rates resulting from fitting analagous spectra as that in (a), 
but over a range of temperatures. 



III. DISCUSSION 

The present data was taken in the 'transverse-field 
geometry', with the field applied perpendicular to the 
initial muon spin directioii^. Previously published 
data on both underdopecP^ and overdoped^^ systems re- 
veal zero relaxation within detection limits when in the 
'longitudinal- field geometry', where the field is parallel 
to the initial muon spin direction. As discussed in detail 
by Savici^^, these longitudinal- field /iSR experiments put 
upper limits on possible transverse fluctuation rates as- 
sociated with the observed magnetism. This observation 
allows us to confidently associate relaxation of the TF- 
/iSR asymmetry with a broadening of the static internal 
magnetic field distribution, which has been an implicit 
assumption in the discussion thus far. 

It is important to note that this broadening of the 
internal field distribution is largely without an accom- 
panying shift. By comparing the precession frequency 
of muons in LSCO to those in a reference crystal of 
Ca(C03), we observe that the average field of the muon 
site is within 0.5mT of the applied field at all temper- 
atures in this study and over the entire doping range. 
This corresponds to a shift of <0.01% for /ioH=5T. In 
contrast, the distribution of fields in Lai.7Sro.3Cu04 has 
a FWHM of 5mT at T=2K in the same applied field. 
Stated simply, muons at different (crystallographically 
equivalent) sites are observing radically different local 
fields at low temperatures, even though the average prop- 
erties are changing very little. A large field- width with 
small average shift is contrary to what one would expect 
for a uniform sample response or commensurate long- 
range order. In fact, quite independent of model, this 
observation implies extensive heterogeneity in the local 
magnetic ground state. This comes as quite a surprise 
in systems which are often touted as canonical Fermi liq- 
uids. It is only appropriate then to speculate on possible 
interpretations for the present data, and the remainder 
of the paper will discuss the feasibility of several such 
models. 

NQR studies have successfully applied a distribution- 
of-hole model to explain line-broadening in underdoped 
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cuprate^^^. In such a picture, holes are quasi-locahzed 
in the copper-oxide planes, and each frequency in a 
given distribution is associated with a unique hole den- 
sity. Though overdoped systems are expected to be more 
metallic (and thus have more delocalized charge carriers), 
it is not unreasonable to expect that the holes in the cur- 
rent samples are likewise distributed. Indeed, this idea is 
seemingly supported by earlier SQUID characterization, 
which revealed a small (< 5%) superconducting volume 
with a Tc ~ 30K in some overdoped samples^. However, 
the shifts in jaSK measurements of LSCO are orders of 
magnitude smaller than those observed by NQR, and the 
sheer magnitude of the present effect makes it difficult to 
explain by a variation of hole densities alone. For exam- 
ple, muons in Lai.7Sro.3Cu04 observe fields on average 
5mT removed from the applied field at low temperature, 
with many muons seeing even larger fields. Yet, there is 
no doping of LSCO which would correspond to a 5mT 
shift. Said another way, there is simply not enough vari- 
ation in shift across the entire LSCO phase diagram to 
explain the distribution of the fields in Lai.7Sro.3Cu04. 
More fundamentally, this model is incomplete, as it does 
not explain the doping evolution of magnetic proper- 
ties nor the source of the doping heterogeneity; it sim- 
ply notes the existence of heterogeneity. Thus, for a 
deeper understanding of these systems, further discus- 
sion of magnetic response is required. 

As mentioned in Sec. I, a similar response to ap- 
plied field was observed by Sonier et al. in crystals of 
LSCO and YBCO, and successfully modeled as a conse- 
quence of heterogeneous superconductivity persisting to 
high temperatures^^. This model was predicated on the 
observation that the rate of relaxation in the underdoped 
YBCO systems was proportional to superconducting Tc, 
and indeed the collective data in underdoped LSCO ap- 
pear to follow a similar trend. However, the field response 
for the three most heavily doped samples increases mono- 
tonically with Sr concentration, despite a weakening of 
superconductivity over the same doping range. Thus, 
superconductivity seems like a less suitable candidate to 
explain a heterogeneous field-response in overdoped sys- 
tems, and implies at least the existence of a separate re- 
laxation mechanism which dominates at for large values 
of X. 

Other /iSR data ha ve be en interpreted as a signa- 
ture of spin magnetism^^ Again, heavily overdoped 
cuprates are known to be free from magnetic order or 
glassiness and are widely expected to be non-magnetic. 
Yet, the Curie term in bulk magnetization data is sug- 
gestive of emergent lo cal momen ts on the overdoped side 
of the phase diagranP^E^I^^I^, and there are theoreti- 
cal grounds to expect such behaviou j^^ l ^-^i It is com- 
monly accepted that the most favourable band for doped 
holes in underdoped systems is the so-called Zhang-Rice 
bancP^ a singlet band formed from the hybridization of 
Cu 3dx2_y2 with O 2px and O 2py orbitals. However, 
long-standing studies by Eskes and Sawatzk};^^ predict 
a triplet band of mostly Cu 3d3z2_j.2 character at only 



slightly higher energies. This triplet band is predicted to 
become more stable with decreasing apical oxygen dis- 
tance, which can be affected by the presence of dopant 
cations^^*^. There are existing x-ray studies that sup- 
port both a decreasing apical oxygen distance^^ ^^^^and 
increased occupation of the Cu 3d322_j.2 orbital^^-^ ^^ l ^^ * in 
overdoped systems. 

The slowing of longitudinal fluctuations amongst such 
emergent moments with applied field could certainly 
broaden the local magnetic field distribution in a sim- 
ilar manner to what is seen. To test the ability of this 
picture to explain our data, we simulated the expected 
/iSR asymmetry function from dipole coupling to a sys- 
tem 0.5/iB aligned moments on 8% of the planar copper 
sites. The density and size of the magnetic moments were 
entirely constrained by bulk magnetization data of Waki- 
moto et a/.E^for Lai.7Sro.3Cu04. As expected for a dilute 
moment system^^, the simulated asymmetry decayed ex- 
ponentially with time. However, even for completely sat- 
urated moments, the maximum expected relaxation rate 
in such a picture is ~1.5/i5~^- about 50% smaller than 
the largest rate seen in Lai.7Sro.3Cu04. More fundamen- 
tally, simulations involving isolated and unconstrained 
(also from bulk magnetization) moments always result 
in a completely isotropic distribution of fields, in stark 
contrast to the large signal anisotropy demonstrated in 
Fig. |4] One may be able to reconcile an isolated mo- 
ment picture with the data, but only if one includes a(n) 
(anisotropic) staggered sample response which at least 
partially compensates a larger local moment. 

Alternatively, Kopp et alW^ have argued that the 
overdoped cuprates may be characterized by competing 
ferromagnetism, with an eventual transition to a fer- 
romagnetic ground state at some critical doping. In 
support of this idea, band calculations by Barbiellini 
and Jarlborg^^ predict that regions of weak ferromag- 
netism should emerge around clusters of dopant ions in 
La2-xBaxCu04 at high dilution. As with the isolated 
moment picture however, it is difficult to reproduce the 
largest observed relaxation rates within a ferromagnetic 
cluster picture, since now the total moment in the sys- 
tem is constrained by bulk magnetization to be 0.5/iB per 
overdoped hole. Specific modeling suggests the ferromag- 
netic cluster model can increase the predicted relaxation 
rate by up to 20% over the isolated moment picture, but 
one is still unable to recreate the size or anisotropy of 
the field-response observed in Lai.7Sro.3Cu04 without in- 
cluding a staggered component to the local magnetism. 

Perhaps a simpler explanation of the observed line- 
broadening is that it reflects clusters of staggered 
magnetic moments, possibly seeded by the dopant 
cations themselves. The temperature, field and 
orientational dependence of the current effect bear 
striking resemblan ce to data previous ly re ported by 

j^Mj^20 2m 23 24 25l26 | 27 | 28 | 29 | 30 J ^SR™ StudicS of 

Zn-doped YBCO, from which many infer the existence 
of staggered islands. Isolated islands with a net mag- 
netic moment would act much like isolated paramag- 
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netic moments, and would broaden the local field dis- 
tribution in similar way. Yet, the net moment of an is- 
land would be greatly reduced from the total local mo- 
ment and could easily reconcile the small effect seen 
in bulk magnetization with the larger effect seen here 
(c.f. Ref. [72]). Specific simulations by us confirm that 
the current data for x=0.30 are consistent with an is- 
land with radius r < one lattice parameter. In such a 
scenario, the linear increase of C{x) with x in Fig. [s] 
is immediately understood to be a consequence of an 
increasing number of nucleation centers for impurity- 
induced magnetism. The rapid decrease in Q{x) upon 
overdoping might reflect a shrinking magnetic correla- 
tion length (and thus island size), which would increase 
the total uncompensated spin of the islands and en- 
hance q=0 magnetic susceptibility at the expense of the 
staggered susceptibility. Consistent with shrinking mag- 
netic correlations, we again mention the disappearance of 
glassy static magnetism^^ and the emergence of Curie- 
like magnetizatioiP^ ^^ * ^^ * ^^ * around the same doping. 
The enhanced field-response of Lai.75Euo.iSro.i5Cu04 
over LSCd^and Yo.8Cao.2Ba2Cu306 over YBCO^^ fur- 
ther supports the central role of dopant cations rather 
than doped holes. 

The interpretation of NMR results in Zn-doped YBCO 
is still being actively discussed^^. However, local re- 
gions of staggered magnetization about dopant Zn^+ 
ions are predicted by a number of theoretical studies 
and expected to arise from correlation effects, even in 
metallic cuprates^^. Andersen et al. recently argued 
that completely analogous physics is expected in Zn- 
free LSCO due to weak Coulombic disorder from Sr^+ 
cation^. Thus, it is perhaps not surprising that we 
have observed data consistent with staggered moments. 
Whether it is appropriate to apply these underdoped 
studies to the more Fermi- liquid- like overdoped systems 
is not obvious and deserves further study. In support 
of this idea however, we note that strong correlations 
have already been inferred from an enhanced Kadowaki- 
Woods ratio in Lai.7Sro.3Cu04^, and both inelastic neu- 
tron scattering^^ and /iSFP^ report an enhancement of 
antiferromagnetic fluctuations with Zn-doping in heavily 
overdoped systems. Antiferromagnetic fluctuations seem 
to decrease with Sr-doping in the Zn-free (and Zn-doped) 
system^^, but this may actually reflect a decreasing mag- 
netic correlation length, inferred above to be less than 
one lattice parameter at x=0.30. Further, the application 
of magnetic field may encourage staggered magnetism 
where little is seen in zero field. Indeed, the central role 
of field is noted in theoretical studies of staggered mag- 
netism in Zn-doped LSCO^, and in at least one cuprate 
system, applied field is seen to modify neutron spectra 
in a way that parallels disordei^^ . To explore these ideas 



more thoroughly, we strongly encourage future measure- 
ments of the inelastic neutron spectrum of heavily over- 
doped cuprate systems in applied field and with varying 
levels of disorder. 

IV. SUMMARY AND CONCLUSIONS 

We have presented TF-/iSR data on overdoped 
La2-xSrxCu04 which shows a large broadening of the 
local magnetic field distribution in response to applied 
field over a range of temperatures. The temperature, 
field and orientational dependence closely reflects similar 
data reported in underdoped materials. Characterization 
of the data through fits to a Curie- Weiss function shows a 
roughly linear increase of C{x) with Sr^+ concentration, 
but a marked shift in Q{x) above optimal doping, where 
a Curie term in the bulk magnetization also emerges. 

Though emergent paramagnetic moments and ferro- 
magnetic clusters have been explored as possible expla- 
nations for the current data, it is difficult to reconcile 
the large local response with the weak bulk magnetiza- 
tion without including a staggered component to the lo- 
cal magnetism. We suggest that a simpler explanation 
of our data is dopant-induced staggered magnetization 
seeded by Sr^+ disorder and stabilized by field. Future 
experiments to investigate the effects of field and disorder 
in the overdoped cuprates are encouraged. 

We emphasize the role that heterogeneity must play in 
any explanation of the current data. This is in contrast 
to the homogeneous Fermi-liquid ground state envisioned 
by many to exist in these systems, and implies at least 
strong correlation effects in the overdoped region of the 
phase diagram. The heterogeneous field response seen 
here complements previous reports of phase separation 
at low temperatures, and suggests that heterogeneous su- 
perconductivity at low temperatures may have its origin 
in a heterogeneous metallic state above Tc. 



V. ACKNOWLEDGEMENTS 

The authors would like to acknowledge useful conversa- 
tions with J.E. Sonier and G.A. Sawatzsky. Work at Mc- 
Master was supported by NSERC and CIFAR. Work at 
Columbia was supported by NSF under DMR-05-02706 
and DMR-08-06846. The work at Lawrence Berkeley 
Laboratory was supported by the Office of Basic en- 
ergy Sciences, U.S. Dept. of Energy under Contract No. 
DE-AC03-76SF0098. We appreciate the hospitality and 
technical assistance of the TRIUMF Centre for Molecu- 
lar and Materials Science, where these experiments were 
performed. 



Current address: Neutron Scattering Science Division, 37831, USA email: macdougallgj@ornl.gov 

Oak Ridge National Laboratory, Oak Ridge, Tennessee 



7 



^ K. M. Lang et aL, Nature 415, 412 (2002). 
2 J. E. Hoffman et aL, Science 295, 466 (2002). 
^ K. McElroy et al., Phys. Rev. Lett. 94, 197005 (2005). 
^ K. McElroy et al., Science 309, 1049 (2005). 
^ J. Lee et al., J. Phys. Chem. Solids 66, 1370 (2005). 
^ T. Kato et al., Phys. Rev. B 72, 144518 (2005). 
^ K. Yoshimura et al., J. Phys. Soc. Jap. 58, 3057 (1989). 
^ H. Ton et al., J. Phys. Soc. Jap. 61, 1477 (1992). 
^ J. Haase et al., J. Supercond. 13, 723 (2000). 
^° J. Haase et al., J. Supercond. 15, 339 (2002). 

P. M. Singer et al., Phys. Rev. Lett. 88, 047602 (2002). 

P. Singer et al., Phys. Rev. B 72, 014537 (2005). 

G. Williams et al., Phys. Rev. B 71, 014515 (2005). 

H. J. Graffe et al., Phys. Rev. Lett. 96, 017002 (2006). 
R. Ofer et al., Phys. Rev. B 73, 012503 (2006). 

R. J. Birgeneau et al. (2006). 

A. T. Savici et al., Phys. Rev. B 66, 014524 (2002). 

K. M. Kojima et al., Physica B 326, 316 (2003). 

H. E. Mohottala et al.. Nature Physics 5, 377 (2006). 
2° H. Alloul et al., Phys. Rev. Lett. 67, 3140 (1991). 
2^ R. E. Walstedt et al., Phys. Rev. B 48, 10646 (1993). 

22 A. V. Mahajan et al., Phys. Rev. Lett. 72, 3100 (1994). 

23 G. V. M. Williams et al., Phys. Rev. B 52, R7034 (1995). 

24 G. V. M. Wilhams and J. L. Tallon, Phys. Rev. B 57, 
10984 (1998). 

2^ J. Bobroff et al., Phys. Rev. Lett. 83, 4381 (1999). 
26 M.-H. Julien et al., Phys. Rev. Lett. 84, 3422 (2000). 
2^ G. Williams et al., Phys. Rev. B 64, 104506 (2001). 

28 Y. Itoh et al., Phys. Rev. B 67, 064516 (2003). 

29 S. Ouazi et al., Phys. Rev. B 70, 104515 (2004). 

^° S. Ouazi et al., Phys. Rev. Lett. 96, 127005 (2006). 
3^ B. Nachumi et al., Phys. Rev. Lett. 77, 5421 (1996). 

32 S. H. Pan et al.. Nature 403, 746 (2000). 

33 G. Williams et al., Phys. Rev. B 71, 176502 (2005). 

34 C. Bernhard et al., Phys. Rev. Lett. 80, 205 (1998). 

3^ B.M. Andersen, P.J. Hirschfeld and M. Schmid, Phys. Rev. 

Lett. 99, 147002 (2007). 
^6 C. Proust et al., Phys. Rev. Lett. 89, 147003 (2002). 
3^ S. Nakamae et al., Phys. Rev. B 68, 100502(R) (2003). 

38 A. Kaminski et al., Phys. Rev. Lett. 90, 207003 (2003). 

39 M. Plate et al., Phys. Rev. Lett. 95, 077001 (2005). 
4° N. Hussey et al., Nature 425, 814 (2003). 

41 C. Panagopoulos et al., Phys. Rev. B 66, 064501 (2002). 

42 Y. Kubo et al., Phys. Rev. B 43, 7875 (1991). 



43 M. Oda et al., Physica C 183, 234 (1991). 

44 T. Nakano et al., Phys. Rev. B 49, 16000 (1994). 

4^ S. Wakimoto et al., Phys. Rev. B 72, 064521 (2005). 
46 H. Mashima et al., Phys. Rev. B 73, 060502 (R) (2006). 
4^ Y. Wang et al., Phys. Rev. B 76, 064512 (2007). 

48 Y. Uemura et al.. Nature 364, 605 (1993). 

49 J. Schutzmann, S. Tajima, S. Miyamoto, and S. Tanaka, 
Phys. Rev. Lett. 73, 174 (1994). 

^° S. Uchida, K. Tamasaku, and S. Tajima, Phys. Rev. B 53, 
14558 (1996). 

^1 H.H. Wen, X. Chen, W. Yang, and Z. Zhao, Phys. Rev. 
Lett. 85, 2805 (2000). 

^2 S. Wakimoto et al., Phys. Rev. Lett. 92, 217004 (2004). 

^3 A. T. Savici et al., Phys. Rev. Lett. 95, 157001 (2005). 

^4 K. Ishida et al., J. Magn. Magn. Mater. 310, 526 (2007). 
J. Sonier et al., Phys. Rev. Lett. 101, 117001 (2008). 

^6 G. J. MacDougall et al., Physica B 374-375, 211 (2006). 
S. Lee, S. Kilcoyne, and R. Cywinksi, eds. (Scottish Uni- 
versities Summer School in Physics and Institute of Physics 
Publishing, Bristol and Philadelphia, 1998). 

^8 G.J. MacDougah et al., Phys. Rev. Lett. 101, 017001 
(2008). 

^9 E.G. Zhang and T. Rice, Phys. Rev. B 37, 3759 (1988). 
^° H. Eskes and G.A. Sawatzky, Phys. Rev. Lett. 61, 1415 
(1988). 

H. Eskes and G.A. Sawatzky, Phys. Rev. B 44, 9656 (1991). 
^2 A. Maignan et al., Physica C 170, 350 (1990). 
^3 C. Martin et al., Chem. Mater. 7, 1414 (1995). 
^4 D. Haskel et al., Phys. Rev. B 56, R521 (1997). 

C. Chen et al., Phys. Rev. Lett. 68, 2543 (1992). 

E. Pellegrin et al., Phys. Rev. B 47, 3354 (1993). 

P. Srivastava et al., Phys. Rev. B 54, 693 (1996). 
^8 A. Kopp, A. Ghosal, and S. Chakravarty, Proc. Natl. Acad. 

Sci. USA 104, 6123 (2007). 
^9 B. Barbiellini and T. Jarlborg, Phys. Rev. Lett. 101, 

157002 (2008). 

P. Mendels et al., Phys. Rev. B 49, 10035 (1994). 

C. Bernhard et al., Phys. Rev. B 58, R8937 (1998). 
^2 p. Mendels et al., Europhys. Lett. 46, 678 (1999). 
^3 M. Gabay et al., Phys. Rev. B 77, 165110 (2008). 
'^4 Risdiana et al., Phys. Rev. B 77, 054516 (2008). 

Y. S. Lee et al., Phys. Rev. B 69, 020502(R) (2004). 



